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SECTION  1 
INTRODUCTION 

A  high  altitude  nuclear  explosion  like  the  1962  burst  Starfish 
(400  km  altitude)  generates  large  currents  in  the  ionosphere.  These  cur¬ 
rents,  like  the  current  in  a  normal  wire,  produce  magnetic  fields  at  points 
far  away  from  the  region  where  the  actual  currents  flow.  When  these  magnetic 
fields  reach  the  surface  of  the  earth  they  set  up  electric  fields  and  cur¬ 
rents  in  the  soil  (or  ocean) .  A  long  conducting  wire  in  this  region  may 
experience  a  large  voltage  drop  across  widely  separated  ground  connections, 
which  can  interrupt  the  line's  function. 

The  source  of  the  ionospheric  currents  is  the  magneto -hydrodynamic 
(MHD)  motion  of  the  ionized  air  surrounding  the  burst  point.  The  term  MHD/EMP 
has  been  coined  to  distinguish  this  effect  from  the  EMP  generated  by  other 
mechanisms.  To  calculate  the  MHD/EMP  signal  one  first  needs  to  calculate 
the  current  patterns  generated  by  a  nuclear  explosion  in  the  ionosphere. 

One  also  needs  to  characterize  the  soil  where  the  long  lines  are  located 
to  relate  the  incident  magnetic  field  to  an  induced  electric  field. 

For  the  past  few  years  the  DNA  MHD/EMP  program  has  been  directed 
at  attacking  these  two  problems.  MRC  has  been  primarily  involved  with  the 
calculation  of  the  current  patterns.  Our  efforts  have  involved  both  trying 
to  develop  simple  models  for  the  ionosphere  currents1  and  actually  trying 
to  compute  the  currents  for  nuclear  events  such  as  Starfish2’3.  The 
modeling  only  met  with  limited  success  because  of  the  complex  nature  of 
the  current  patterns.  However,  it  did  lead  us  to  conclude  that  there  were 
several  distinct  mechanisms  responsible  for  the  signals  measured  after  the 
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high  altitude  events  during  the  1962  test  series  (see  the  Appendix).  In 
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particular,  Longmire4  concluded  that  the  700  gamma  (7x10"  gauss)  signal 
observed  about  70  seconds  after  the  Starfish  event  was  due  to  the  upward 
motion  (heave)  of  the  atmosphere  below  the  burst  point.  This  air,  which  is 
initially  between  90  and  150  km  altitude,  is  heated  by  x-rays  and  by 
energetic  ions  which  travel  down  the  magnetic  field  lines  from  the  burst 
point.  These  regions  of  heated  air  are  called  the  x-ray  patch  and  the 
kinetic  energy  (or  debris)  patch. 

The  motion  of  the  air  after  the  burst  is  computed  using  the  3-D 
MHD  code  MICE5.  MICE  is  designed  primarily  as  a  radar  and  communication 
effects  code  so  the  emphasis  is  on  correctly  computing  electron  densities. 
Electric  and  magnetic  fields  are  computed  since  they  influence  the  motion 
of  the  plasma,  but  with  approximations  which  produce  somewhat  inaccurate 
results  over  much  of  the  grid.  Previous  attempts  to  use  MICE  to  compute 
MHD/EMP  effects  were  generally  unsuccessful2’3.  For  the  present  effort 
we  spent  some  time  trying  to  change  MICE  to  do  a  better  job  at  computing 
the  magnetic  fields  but  finally  decided  on  a  new  approach. 

To  solve  for  the  MHD/EMP  at  the  surface  of  the  earth  we  use  the 
neutral  wind  velocities  and  the  ion  and  neutral  densities  computed  from 
MICE  (quantities  it  was  designed  to  compute  accurately)  as  input  to  a 
SIMPLE  3-D  Maxwell  equation  solver.  The  Maxwell  solver  (MHDEMP)  is  optimized 
to  solve  for  the  electric  and  magnetic  fields  that  result  from  the  atmospheric 
motion  computed  by  MICE. 

Using  this  combination  of  MICE  and  the  MHDEMP  code  we  were  able, 
for  the  first  time,  to  compute  accurate  magnetic  fields  at  the  surface  of 
the  earth  for  the  Starfish  event.  We  obtained  excellent  agreement  with 
the  measured  signal  at  Johnston  Island  (see  Figure  12,  pg  27).  It  is 
dangerous  to  conclude  from  one  measurement  that  our  analysis  is  correct 
at  all  other  locations  but  before  the  technique  was  developed  there  was 
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not  even  agreement  at  the  one  point.  Wt  have  demonstrated  that  we  understand 
the  mechanism  responsible  for  the  large,  late  time,  magnetic  pulse  seen 
after  Starfish  and  that  we  can  compute  the  observed  signal  from  first 
principles.  As  we  have  time  to  study  the  computer  output  and  increase  our 
understanding  of  the  phenomena  we  will  be  able  to  develop  more  confidence 
in  the  results  at  other  points  and  for  other  bursts. 

In  Section  2  we  briefly  explain  the  operation  of  the  MHDEMP  code. 
Since  this  section  will  only  be  of  interest  to  someone  familiar  with  the 
high  altitude  MHD  equations  we  have  assumed  a  working  knowledge  of  the 
terms  and  equations  generally  used  in  that  field.  (See  Chapter  12  of  Ref¬ 
erence  6  for  example.) 

Three  separate  MHD/EMP  simulations  were  made  for  this  report.  The 
first  simulation  used  the  standard  MICE  plasma  parameters  as  input.  This 
gave  the  best  results  and  is  presented  in  Section  3.  The  basic  MICE  calcu¬ 
lation  does  not  include  the  ionization  produced  by  g-particles  which  are 
emitted  by  the  bomb  debris  and  spiral  down  the  magnetic  field  lines  to  an 
altitude  of  60  to  80  km  making  what  is  called  a  beta  patch.  It  was  felt 
that  this  ionization  should  be  important  to  the  MHD/EMP  signal  so  an  addi¬ 
tional  run  was  made  which  included  the  beta  patch.  Also  there  was  some 
thought  that  the  MICE  electron  densities  might  be  as  much  as  a  factor  of 
three  too  high  so  a  third  simulation  was  run  with  the  MICE  conductivities 
artificially  set  a  factor  of  three  smaller.  A  beta  patch  was  also  included 
for  this  third  simulation.  The  results  of  these  two  additional  simulations 
are  included  in  Section  4. 

A  summary  of  the  results  is  presented  in  Section  5  along  with 
recommendations  for  future  work.  We  have  also  included  in  the  Appendix  a 
brief  time  history  of  events  following  a  high  altitude  nuclear  burst 
relevant  to  the  MHD/EMP  problem.  A  reader  not  familiar  with  all  the 
phenomena  involved  might  benefit  from  reading  the  Appendix  before  reading 
the  following  sections. 
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SECTION  2 

THE  MHDEMP  CODE 


Although  it  does  not  appear  feasible  to  use  the  MICE  currents 
directly  to  compute  B  on  the  ground,  we  believe  the  other  plasma  parameters 
from  MICE  are  very  reliable.  The  new  MHDEMP  code  computes  B  on  the  ground 
using  the  neutral  heave  and  conductivities  from  MICE.  The  MHDEMP  code  is  a 
straightforward  3-D  code,  which  can  be  rapidly  iterated  many  times.  The 
grid  extends  to  the  ground  so  the  B-field  there  can  be  computed  directly. 

Our  immediate  objective  was  to  use  the  co-e  to  calculate  the  large  negative 
signal  measured  on  the  ground  70  seconds  after  the  Starfish  event.  It 
appears  the  code  is  capable  of  computing  the  earlier  positive  signal,  as 
well . 


TB»  MHDEMP  code's  basic  assumption  is  that  the  electromagnetic 
forces  balance  the  neutral  drag  force.  This  is  the  same  basic  assumption 
used  for  the  MICE  code  and  other  MHD  codes,  except  that  usually  other  forces 
are  included  in  addition  to  the  neutral  drag  force.  The  MICE  code  includes  iner 
tia,  gravity,  and  the  pressure  gradient  force.  These  forces  should  not  be  impor 
tant  for  our  times  of  interest.  The  resultant  equation  for  the  current  is6 


J  =  o  E  +  a 

ii  ii  p 


(e  ) 

\  x  c  /  H\  c  nx/ 


(1) 


where  we  have  introduced  the  three  well  known  conductivities: 


eN 

Parallel  a  = -g- 


l  +  n  /  n . 

eo  10 


n  +  n  . (i  +  n  /n.  )  , 

eo  ei  eo  10  J 


(2) 
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Pedersen  a 


(n-  +n  ) [i+n-  n  +  q  .  (q.  +q  )1 

10  eo  L  10  eo  ei  10  eo  J 

2  2 

(l+n-  ) (i+n  +  n  .  (q.  +n  )f2+2n.  n  +n  .(n.  +n  )] 

io  eo  ei  10  eo  10  eo  ei  *■  10  eo^J 


Hall  aH  =  ^  - j - 2 - ~ ^ - 

Ld+nio)(i+neo)  -  nei(nio+neo)[2+2nioneo+nei(nio+neo)]. 


The  n's  are  dimensionless  ratios  of  gyro  frequency  to  momentum  transfer 

collision  frequency:  q.  =  v.  /Q. ,  q  =  v  /Q  and  q  ■  =  V  . /Q  . 

io  io  i’  ‘eo  eo  e  'ei  ei  e 


One  can  invert  Equation  1  to  obtain  E. 


E  =  -L  J 

"  a„  " 


,  /  _  9E  \ 

^((VX%  -  f  Iff") 


Ex  =  2  2 

(V°H} 


r  2  2. 

(W 


/  ^  „  9E  \  V  > 

,  ((,XBh  -  C  St)  -  -T 


'4,(0 

pH 

The  second  lines  of  Equations^  and  6  were  written  using  one  of  Maxwell's 
Equations,  J  =  (VxB  -  2-  ^  ) .  The  second  line  of  Equation  6  ig¬ 
nores  the  symmetry  breaking  Hall  term,  which  reflects  the  present  state 
of  the  MHDEMP  code.  This  term  is  believed  to  be  an  insignificant  factor 
for  most  of  the  grid,  but  may  be  important  enough  from  80  km  to  180  km 
altitude  to  warrant  its  inclusion  into  the  code.  Its  omission  allows  the 
use  of  a  symmetry  plane  and  better  resolution  in  the  MHDEMP  code. 
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This  paper  uses  modified  Gaussian  cgs  units.  The  current  density 

-2 

and  magnetic  fields  are  in  emu,  abamps  cm  and  Gauss,  electric  fields 
in  esu,  statvolts  cm  \  The  dimensions  of  conductivity  are  cm  *  and  the 
dielectric  constant  is  dimensionless. 


The  code  uses  Equations  5  and  6  to  find  E 
then  advances  B  in  time  via 


from  B  and  V 

n 


and 


9B 

3t 


cVxE 


(7) 


To  speed  execution  time  the  MHDEMP  code  is  very  straightforward.  For 
#  -+ 

instance,  the  parallel-to-B  direction  for  the  conductivities  is  taken  to 

be  the  ambient  B  direction.  It  would  take  a  large  computational  effort  to 

compute  the  conductivity  tensor  with  the  true  parallel-to-B  direction,  which 

normally  is  close  to  the  ambient  direction.  The  conductivities  and  (V  *B/c) 

n 

terms  are  updated  only  at  MICE  dump  times,  which  are  spaced  at  5  to  15  sec¬ 
ond  intervals.  These  simplifications  may  be  eliminated  in  the  future,  if 
more  accuracy  is  desired. 


The  MHDEMP  code  will  be  unstable  if  the  signal  speed  is  too  fast 
compared  to  the  time  step.  The  fastest  signals  occur  in  the  lower  atmosphere 
from  the  ground  up  to  the  altitude  where  significant  electron  density  occurs, 

which  is  either  the  beta  patch  ( _ 60  km)  or  the  E-region  (~100  km).  The 

conductivity  of  the  lower  atmosphere  is  insignificant  and  the  signal  speed 
reduces  to  the  speed  of  light  in  a  dielectric,  c/v/e.  The  dielectric  constant 
for  the  entire  computational  grid  is  artificially  increased,  so  that 


c(At) 

ASq(.8) 


where  ASq  is  the  smallest  cell  dimension  anywhere  in  the  MHDEMP  grid  and 
At  is  the  time  step.  Typical  operating  parameters  are  ASq  =  8  km. 
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At  =  0.02  seconds  and  e  =  9*10^.  The  early  development  versions  of  the 

£  8  E 

MHDEMP  code  omitted  the  displacement  current  (-—•»—  in  Equation  6  and  7) 

C  dt 

and  obtained  stability  by  introducing  an  artificial  conductivity.  The 
displacement  current  is  normally  omitted  from  MHD  codes.  Without  the 
displacement  current  the  signal  speed  is  the  3-D  magnetic  diffusion  speed 
and  the  stability  requirement  is 


a 

ii 


and 


V°H  j  ^  cAt 
QP  '  (1/6)  4tt  (AS) 2 


(9) 


2  2 
/  °P+aH  \ 

The  code  checked  the  input  'alues  of  a  and  ^ — - — J  in  every  cell  and 
increased  tl  cm  to  the  minimum  stable  value,  if  necessary.  The  AS  is  the 
smallest  of  the  three  cell  dimensions  Ax,  Ay  or  Az.  This  method  is  analogous 
to  the  Boris  mass  modifications 7  of  the  MICE  code  (where  the  signal  speed  is 


the  Alfven  speed)  and  also  results  in  artificial i>  high  currents  in  regions 
where  the  conductivities  should  be  low. 


The  artificial  conductivities  proved  to  be  a  significant  liability. 
Not  only  were  the  currents  artificially  increased  in  peripheral  regions  of 
the  grid,  but  a  sizeable  amount  of  current  was  able  to  flow  through  the 
lower  atmosphere  into  the  bottom  boundary  of  the  grid.  This  introduced 
spurious  magnetic  fields  on  the  ground  comparable  in  magnitude  to  the 
expected  fields.  Decreasing  the  time  step  and  adding  a  special  resistive 
layer  helped  slightly,  while  increasing  computer  costs  a  great  deal.  For 
this  reason  we  introduced  the  displacement  current.  The  conduction  current 
flowing  through  the  lower  atmosphere  is  now  a  factor  of  10^  less  than  in 
the  early  versions  of  the  code. 

The  location  of  the  top  boundary  is  chosen  so  that  no  signal 
generated  above  the  top  boundary  has  time  to  diffuse  to  the  ground  by 
the  end  of  the  simulation.  Presently,  the  perfectly  conducting  ground 
and  the  symmetry  plane  use  the  boundary  condition 
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(10) 


n  x  J  =  0  . 


The  top  and  other  side  boundaries  obey 
A 

n  •  J  =  0 

to  keep  current  from  leaving  the  grid. 


(11) 


The  beta  patch  is  ignored  by  MICE,  since  it  lies  below  the  bottom 
of  the  MICE  grid.  The  MHDEMP  code  uses  a  beta  patch  model  previously  used 
by  Crevier  combined  with  some  features  of  a  model  by  Longmire8.  The  power 
into  betas  is 


(.15)  ( . 05 ) oof 

p  =  - - —  ergs/sec  (12) 


where  is  the  fission  yield,  the  .05  represents  the  5%  of  fission  yield 

into  betas  and  the  0.15  is  a  normalizing  factor.  The  rate  of  ion-pair 
production  in  the  beta  patch  is 


n 


(1.7  x  35  x 1.6X10'12} 


-3 

cm 


-1 

sec 


(13) 


The  assumptions  that  lead  to  Equation  13  are  as  follows.  The  beta  energy 
is  taken  to  be  1.7  MeV.  The  energy  loss  rate  is  (1.6  MeV  cm2/gram) p^. 

An  ion-electron  pair  is  formed  for  every  35  eV  deposited  by  the  betas. 

One  half  of  the  betas  go  up  and  one  half  go  down.  The  /2  comes  from 
the  assumption  that  the  betas  spiral  around  the  field  lines  with  a  pitch 
angle  of  45°.  The  large  pitch  angle  is  due  to  the  magnetic  mirror  effect, 
since  the  betas  travel  from  a  region  where  the  fields  are  weak  to  the 
ambient  field. 
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a 


* 


i 

i 


fg(cm“2) . 


The  relative  number  of  betas  on  each  field  line  is  determined  by 
Longmire's  model8  for  f^  when  applied  to  Starfish  yields 


f  =  t/l  +  1 . 64 (S/ S  )  e"S/S°  (14) 

6  4ttS2  0 

o 

where  S  is  the  distance  from  the  burst  point  field  line  to  the  field  line 
in  question  and  Sq  is  60  km. 


The  equilibrium  electron  density  is  determined  from  n  by  lumping 
the  deionization  chemistry  into  a  single  recombination  rate 


n  =  ;  a  =  3*10~7  cm3/sec  (15) 

The  peak  ionization  occurs  at  60  km  where  the  bulk  of  the  beta  energy  is 
deposited.  Ionization  due  to  betas  above  100  km  or  below  60  km  is  ignored. 
To  compute  conductivities  below  the  bottom  of  the  MICE  grid  (77  km),  ambient 
electron  temperatures  and  neutral  densities  are  assumed.  The  grid  spacing 
for  the  Starfish  run  described  in  Section  3  was  8  km  in  the  beta  patch 
region,  which  is  too  coarse  to  resolve  the  beta  patch  at  the  ionization 
peak.  A  finer  grid  spacing  would  be  better. 
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SECTION  3 

STARFISH  SIMULATION 


This  section  describes  a  Starfish  simulation  using  the  MHDEMP  code. 
Starfish  was  chosen,  since  among  the  Fishbowl  events,  it  had  the  largest  sig¬ 
nal  and  represents  the  greatest  threat  to  long-line  systems.  In  addition,  a 
recent  MICE  Starfish  calculation  had  just  been  completed  out  to  105  seconds. 

F.  E.  Fajen  kindly  made  these  results  available  to  us.  While  we  were  in  the 
midst  of  performing  the  MHDEMP  simulations,  a  new  MICE  Starfish  run  was 
announced  for  later  this  year.  It  is  anticipated  that  the  new  Starfish 
results  could  have  a  factor  of  3  lower  electron  densities  in  the  debris 
patch.  Our  first  MHDEMP  simulation  used  the  existing  MICE  plasma  parameters 
and  the  beta  patch  model  explained  in  Section  2.  The  second  simulation  had 
the  same  beta  patch  and  MICE  parameters  except  the  MICE  electron  densities 
were  decreased  by  a  factor  of  3.  This  was  done  to  check  the  sensitivity  of 
the  results  to  Ng  in  anticipation  of  the  future  MICE  Starfish  run.  The 
third  simulation  was  the  same  as  the  first,  except  the  beta  patch  was  excluded. 
The  third  simulation  gives  the  best  fit  to  the  J.I.  data  and  will  be  described 
in  depth.  The  results  of  the  first  two  simulations  will  be  discussed  in 
Section  4.  All  figures  in  this  section  will  come  from  the  third  simulation, 
i.e.,  without  the  beta  patch. 

The  simulations  starts  at  20  seconds  with  B  ambient  everywhere. 

The  MHDEMP  grid  is  the  same  as  the  MICE  grid  (+  x  is  north;  +  y  is  west, 

+z  is  up),  except  the  bottom  of  the  MHDEMP  is  at  z  =  0  rather  than 
z  =  77  km  and  the  top  is  at  z  =  400  km  rather  than  z  =  1024  km.  Signals 
generated  above  400  km  are  not  expected  to  propagate  to  the  ground  by  the 
end  of  the  simulation  at  110  seconds.  The  origin  of  the  coordinate  system 
(x=0,  y=0,  z=0)  is  located  on  the  ground  directly  beneath  the  burst  point. 
•Johnston  Island  is  at  (x  =  30  km,  y  =  0,  z  =  0) .  The  time  step  was 
0.02  seconds  and  e  =  9x10^  was  needed  for  stability. 
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Page  44  shows  Longmire's  prediction  of  the  current  patterns  in 
the  ground  and  in  the  debris  patch,  which  agrees  for  the  most  part  with  the 
result  of  the  MHDHMP  code.  The  rising  neutrals  operate  as  a  homopolar  gener¬ 
ator.  From  Equation  1  one  would  expect  a  westward  current  in  response  to 

the  V  xB/c  term.  This  current  must  find  an  eastward  return  path.  The 
n 

current  will  diminish  or  stop  depending  upon  the  magnitude  of  an  opposing 

E  ,  if  such  an  E  is  able  to  be  established.  The  current  is  expected  to 

1  ±  2  2 
flow  along  field  lines  in  regions  where  a  »  ((a  +a  )/o  )  and  is  able  to 

2  9  11  ^  ^  * 

cross  field  lines  easily  where  (0p+a^)/ap  is  large.  Therefore,  the  illus¬ 
tration  on  page  44  shows  the  current  flowing  predominately  down  field  lines 

2  2 

and  returning  to  the  east  at  lower  altitudes  where  (0^+0^) /op  is  large. 

The  illustration  shows  another  loop  to  the  south  which  is  above  the  burst 
point  and  originally  was  meant  to  represent  the  southern  conjugate  point. 

Our  simulation  had  a  similar  current  pattern  in  the  debris  patch  and  on  the 
ground  with  some  exceptions.  The  southern  current  loop  was  not  above  the 
generator  region,  but  was  below  the  generator  region.  Also  there  exists 
another  load  region  not  predicted  in  the  illustration  on  page  44. 


Figure  1  shows  a  sketch  of  the  perturbed  field  AB  on  the  ground 
at  69.2  seconds.  Also  shown  is  the  direction  of  the  surface  currents  induced 
in  the  ground  to  keep  the  perturbed  magnetic  field  from  penetrating  the  per¬ 
fectly  conducting  ground.  The  surface  current  is  perpendicular  to  AB  and 
flows  in  two  opposing  rings.  For  an  earth  with  a  large  finite  conductivity 
the  values  of  the  surface  current  would  be  computed  via  Stokes's  Law  from 

| AB | ,  if  the  skin  depth  is  known.  Letting  AB  =  0  one  skin  depth  inside 
— ► 

the  earth  and  if  AB  is  parallel  to  the  earth's  surface. 


0> 


jAB[ 

4TTh 


(16) 


where  {j)  is  the  current  density  averaged  over  a  skin  depth  h  and 
| AB  j  is  evaluated  at  the  earth's  surface.  Figures  2  to  4  shows  | AB 
on  the  ground  at  various  times.  We  expect  the  difference  between  a 


perfectly  conducting  boundary  and  a  realistic  conducting  ocean  should  not 
affect  these  values  appreciably. 

The  magnetometer  at  J.I.  measured  the  change  in  the  magnitude  of 
the  total  magnetic  field  B,  ambient  plus  perturbed.  Figures  5  to  7 
show  the  change  relative  to  ambient  of  the  total  field  on  the  ground 
A|"5|  =  +A$|  -  !  B  | .  The  early  times  are  influenced  by  the  fact  that 

the  simulation  started  at  20  seconds  with  ambient  B^.  The  pattern  of  the 
field  on  the  ground  sketched  in  Figure  1  sets  in  by  40  seconds  and  remains 
unchanged  for  the  rest  of  the  run.  In  the  central  region  where  AB  points 
generally  south,  |B|  is  decreased.  On  the  south  (x  =  -  600  km  to  x  =  0) 
and  the  north  (x  =  600  km  to  x  =  1000  km)  ends  of  the  grid,  AB  points 
generally  north  and  |b|  is  increased. 

Between  the  regions  of  positive  A|b|  and  negative  A|B|  are 
two  node  lines.  Magnetometers  sitting  on  a  node  would  indicate  no  change 
in  the  magnitude  of  the  magnetic  fields.  J.I.  is  close  to  the  southern 
node  line. 


From  the  prediction  on  page  44  one  would  expect  to  see  a  two 
ring  current  pattern  in  the  debris  patch.  The  actual  3-D  geometry  is  more 
complicated  than  the  simple  discrete  circuit  pictured  on  page  44.  In  our 
simulation  there  is  no  sharp  distinction  between  the  northern  and  southern 
load  regions.  Figures  8  and  9  show  Jy  in  the  symmetry  plane.  There  is  a 
load  region  (eastward  current)  immediately  below  the  main  generator  region 
due  to  slower  upward  winds  there.  This  load  was  not  predicted  on  page  44. 
There  is  another  load  region  in  the  H-region  at  100  km  altitude. 

Both  the  southern  and  northern  loads  are  at  this  altitude.  Some  of  the 
current  flows  directly  east  into  the  symmetry  plane  after  reaching  the 
F-region.  This  major  return  path  centered  around  x  =  400  km  corresponds 
directly  to  the  northern  return  path  on  page  44.  Since  this  region  is 
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Figure  8.  in  the  symmetry  plane  at  24.6,  40.0,  60.0  seconds. 
MICE  N  without  beta  patch. 
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Figure  9.  3  in  the  symmetry  plane  at  69.2,  90.0,  110.0  seconds 

MICE  N  without  beta  patch. 


a  conducting  plane,  the  current  tends  to  spread  out  as  it  flows  east.  At 
40  seconds  it  is  evident  that  significant  current  is  flowing  into  the 
symmetry  plane  at  a  point  further  south  than  any  of  the  significant  gener¬ 
ator  currents.  By  69.2  seconds  this  level  of  current  extends  to  the  south¬ 
ern  edge  of  the  grid,  and  at  the  same  time  the  magnitude  of  the  northern  load 
currents  have  noticeably  diminished  in  the  region  100  km  <  x  <  700  km. 

The  strength  of  the  southern  return  currents  appears  to  increase 
slightly  during  the  run  compared  to  the  northern  return  currents.  The 
pattern  of  A|b|  on  the  ground  also  shifts  slightly  to  the  south.  Inspection 
of  Figures  5  to  7  shows  this  to  be  the  case.  It  is  not  a  dramatic  change, 
but  the  southern  node  line  shifts  from  x  =  100  km  to  x  =  0  km  from  40  seconds 
to  69.2  seconds.  This  shift  does  have  a  significant  effect  for  magneto¬ 
meters  operating  at  or  around  J.I.  Ten  observer  points  where  chosen  along 
the  symmetry  plane  for  time  plots  of  A|b|  .  Figure  10  is  for  an  observer 
at  x  =  -  184  km.  The  wiggles  in  the  time  plot  are  artificial  electromagnetic 
waves  in  the  lower  atmosphere  bouncing  off  of  the  side  walls  of  the  grid  and 
are  not  real.  We  have  varied  the  conductivities  in  the  lower  atmosphere  and 
the  dielectric  constant  and  verified  that  these  waves  do  not  alter  the 
mean  A|§|.  The  observer  at  x  =  -  184  km  always  sees  a  |B|  greater  than 
ambient  consistent  with  the  fact  that  the  node  line  does  not  drift  this  far 
south.  Figure  11  shows  A  J  B |  for  an  observer  at  x  =  23.4  km,  the  closest 
observer  to  J.I.  The  node  passes  through  this  point  giving  an  early  positive 
signal  followed  by  a  negative  signal.  This  is  the  correct  sequence  of 
events,  but  the  comparison  to  actual  data  is  poor.  Figure  12  shows  the 
next  observer  at  x  =  179.2  km  north.  This  observer  is  always  north  of  the 
node  and  sees  a  negative  signal,  which  is  an  amazingly  good  fit  to  the  J.I. 
data.  The  excellence  of  the  fit  can  not  be  taken  totally  seriously,  since 
our  input  data  is  not  expected  to  have  better  than  a  factor  of  two  type 
accuracy  and  perhaps  is  worse  considering  a  new  Midi  Starfish  is  being 
planned.  Taking  Figures  11  and  12  together  indicates  that  the  MHDEMP 
code  probably  handles  the  essential  details  correctly. 
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|^Bj  (Gammas) 
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We  have  only  discussed  the  two  load  regions  at  100  km  altitude 
and  have  ignored  the  generator  region  and  the  load  region  immediately 
below  it.  During  the  simulation  the  orientation  of  the  generator  and  this 
load  region  changes  and  they  both  rise.  This  should  also  influence  the 
change  in  the  pattern  of  the  field  on  the  ground. 

Our  results  show  that  J.I.  is  at  a  particularly  critical  location. 
The  positive  pulse  followed  by  a  negative  pulse  occurs  only  in  a  limited 
region  on  the  ground.  It  would  require  extremely  accurate  input  data  to 
reproduce  the  correct  magnitude  and  time  of  both  the  positive  and  negative 
pulses  and  have  them  located  precisely  at  J.I. 
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SECTION  4 

ADDITIONAL  CALCULATIONS 

This  section  describes  the  differences  between  the  three  MHDEMP 

simulations.  The  first  simulation  used  the  standard  MICE  plasma  parameters 

as  input  and  had  a  beta  patch.  The  second  simulation  used  a  third  of  the 

MICE  electron  density  (denoted  by  N  /3),  but  had  the  same  beta  patch  as 

6 

the  first  simulation.  The  third  simulation  used  the  MICE  N  ,  but  omitted 

e’ 

the  beta  patch.  The  third  simulation  was  described  in  detail  in  the 
previous  section. 

It  is  not  possible  to  get  a  clear  picture  of  the  large  scale 
differences  between  the  three  runs  from  the  A | B j  observations  around  J.I., 
since  J.I.  is  strongly  influenced  by  the  nearby  presence  of  the  node  line. 
The  first  two  runs  have  nearly  identical  time  histories  of  A | B |  at  J.I. , 
but  are  different  at  all  other  observer  points.  The  field  near  J.I.  appears 
to  change  dramatically  during  the  runs,  while  the  change  to  the  overall  AB 
pattern  is  just  a  slight  shift  to  the  south. 

The  peak  A | B |  on  the  ground  is  a  simple  measure  of  the  important 
large  scale  differences  between  the  runs,  as  well  as  being  an  important 
parameter  for  system  EMP  effects.  The  time,  location  and  magnitude  of  the 
peak  A|?|  on  the  ground  is  shown  in  Table  1.  The  data  in  Table  1  is 
taken  from  the  ten  observer  points  on  the  symmetry  plane.  The  first  two 
runs  have  a  beta  patch  and  therefore  have  more  conductivity  between  the 
kinetic  energy  patch  of  the  nuclear  burst  and  the  ground  than  the  third 
run.  The  third  run  without  a  beta  patch  showed  that  A|B|  peaked  earlier 


Table  1.  Comparison  of  peak  signal  on  ground. 


IH9JUI 

Time 

Location 

1.  Standard  MICE 

e 

Standard  Beta  Patch 

-700y 

110  seconds 

■a 

2.  MICE  Ng/3 

Standard  Beta  Patch 

-580y 

110  seconds 

x  =  300  km 

3.  Standard  MICE  N 

e 

No  Beta  Patch 

-950y 

70  seconds 

■ 

since  the  magnetic  diffusion  time  to  the  ground  was  shorter  due  to  the 
lesser  conductivity.  The  magnitude  of  A[b|  for  the  third  run  was  greater 
since  there  was  less  shielding  by  the  conductivity. 

As  expected,  the  second  run  with  only  one-third  of  the  MICE  Ng 
had  less  current  in  the  generator  region  and  therefore  the  smallest  signal 
on  the  ground.  The  computed  variation  in  the  magnitude  of  the  peak  A | B | 
among  the  three  runs  is  smaller  than  one  might  have  expected. 

The  location  of  the  peak  A | "S |  occurred  farther  to  the  south 
for  the  first  two  runs.  The  beta  patch  increases  the  conductivity  more 
in  the  northern  load  region  than  in  the  southern  load  region.  Figure  13 
shows  the  load  and  generator  currents  for  the  standard  MICE  Ng  run  with 
the  beta  patch.  At  40  and  69.2  seconds  the  return  currents  in  the  con¬ 
ducting  plane  appear  to  be  just  as  spread  out,  if  not  more  so,  as  those 
in  the  run  without  the  beta  patch  shown  in  Figures  8  and  9.  At  110  seconds 
the  return  currents  have  shifted  to  the  north.  This  run  had  the  same  two 
ring  current  patterns  on  the  ground  and  Figure  14  shows  that  the  southern 
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node  line  drifted  south  more  than  in  the  run  without  the  beta  patch.  The 
fact  that  the  node  moved  south  can  not  be  associated  with  the  current  pattern 
in  the  E-layer  shifting  south,  since  they  went  north.  At  110  seconds  the 
southern  current  loop  on  the  ground  is  practically  gone. 

The  Ng/3  run  had  insignificant  current  in  the  southern  load  region 
and  only  had  one  current  ring  on  the  ground  and  only  the  northern  node  line. 
Figure  15  shows  A|b|  on  the  ground  for  the  Ne/3  run. 
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SECTION  5 

CONCLUSIONS  AND  POSSIBLE  FUTURE  IMPROVEMENTS 


Some  questions  are  still  to  be  resolved.  We  only  had  one  current 
loop  for  the  Ng/3  run.  We  are  not  sure  exactly  why  the  A  j  B |  pattern  drifts 
south  during  the  run.  It  is  probably  related  to  the  geometry  of  the  currents 
charging  as  the  generator  region  gains  altitude.  The  model  without  the  beta 
patch  gives  a  better  fit  to  the  J.I.  data  than  the  beta  patch  model.  The 
better  fit  may  be  accidental. 

The  symmetry  breaking  Hall  term  should  be  incorporated  into  the 
MHDEMP  code  on  an  experimental  basis.  It  should  make  a  difference  in  the 
load  region  near  100  km  altitude.  Our  beta  patch  model  may  be  too  rough, 
since  we  are  using  8  km  grid  spacings  in  altitude.  It  may  be  that  we  are 
overestimating  the  beta  patch  or  the  effective  beta  patch  conductivi¬ 
ties  may  be  less  than  imagined  due  to  irregularities  in  electron  density 
caused  by  early  time  instabilities  in  the  debris-air  mixing.  The  beta 
patch  was  observed  to  be  irregular.  The  effective  conductivity  through 
the  beta  patch  should  be  derived  by  averaging  the  resistivities  along 
the  current  path.  This  would  lead  to  a  smaller  effective  conductivity 
than  our  smooth  distribution  of  beta  patch  electrons. 

Forcing  the  currents  to  close  inside  the  MICE  grid  may  affect  the 
fine  features  of  the  AB  pattern  on  the  ground,  such  as  the  precise  loca¬ 
tion  of  the  node  line.  This  is  true  for  the  currents  in  the  E-layer,  since 
this  region  extends  around  the  world.  Currents  in  the  kinetic  energy  patch 
appear  far  enough  away  from  the  boundaries.  If  more  accuracy  is  desired, 
the  MHDEMP  grid  could  be  expanded  beyond  the  MICE  grid.  The  plasma  para- 


37 


meters  are  updated  only  at  MICE  dump  times,  which  are  spaced  15  seconds 
apart  by  the  end  of  the  simulation.  Better  accuracy  could  be  obtained  by 
interpolating  the  parameters  in  time  and  updating  more  often.  This  would 
also  decrease  the  magnitude  of  the  artificial  electromagnetic  waves  in  the 
lower  atmosphere  which  are  responsible  for  the  wiggles  in  Figures  10  to  12. 
These  waves  are  excited  whenever  the  plasma  parameters  are  updated.  Future 
MHDEMP  simulations  should  be  started  at  earlier  times  than  20  seconds  to 
reproduce  the  positive  signal.  In  that  case  it  may  be  worthwhile  to  include 
the  pressure  gradient  as  well  as  the  neutral  slip  as  MICE  does.  These  simu¬ 
lations  were  started  at  20  seconds  since  we  originally  were  concentrating 
on  the  late  time  negative  signal. 

In  summation,  the  MHDEMP  code  working  with  MICE  input  data 
reproduced  the  general  features  of  the  Starfish  data.  IVe  are  now  ready 
to  use  other  existing  MICE  runs  with  various  yields  and  burst  altitudes 
to  determine  which  bursts  produce  the  greatest  magnetohydrodynamic  EMP 
effects  on  long-line  systems. 
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APPENDIX 

THE  MAGNETOHYDRODYNAMIC  EMP 
BRIEFING  PREPARED  FOR  CAPT.  MIKE  BELL,  DNA  (RAEV) 

BY 

CONRAD  L.  LONGMIRE,  MRC 
October  1  978 

1.  HIGH-ALTITUDE  NUCLEAR  EXPLOSIONS  GENERATE  ELECTROMAGNETIC  FIELD  (OR  EMP)  THROUGH 
SEVERAL  PROCESSES. 

1.1  GAMMA  RAYS - ►CURRENTS  OF  COMPTON  RECOIL  ELECTRONS - ►EM  FIELDS. 

GAMMA  INDUCED  EMP  IS  LARGE  FOR  ABOUT  1  ySEC,  THEN  FALLS  OFF.  SEEN  BY 
OBSERVERS  FROM  GROUND  LEVEL  UP. 

THIS  IS  THE  MOST  WELL  KNOWN  MECHANISM. 

THEORY  ESTABLISHED  IN  1964. 

COMPUTER  CODES:  CHAP  AND  HEMP  B  FOR  EARLY  TIMES,  LHAP  AND  HAPS  FOR 
LATE  TIMES. 

CALCULATIONS  ARE  IN  GOOD  AGREEMENT  WITH  EXPERIMENTAL  DATA  FROM  OPERATION 
FISHBOWL. 

1.2  X  RAYS - ►PHOTOELECTRONS - ►EM  FIELDS. 

THIS  EFFECT  IS  PRIMARILY  A  NON-PROPAGATING  PLASMA  OSCILLATION,  SEEN  IN 
THE  REGION  ABOVE  50-KM  ALTITUDE  WHERE  X  RAYS  ARE  ABSORBED,  BELIEVED 
INSIGNIFICANT  FOR  LOWER  OBSERVERS. 

SOME  THEORY.  CODES:  XEMP,  PEMP. 

1.3  NEUTRONS - ►COLLISIONS  WITH  AIR  ATOMS - ►GAMMA  RAYS - ►EM  FIELDS. 

SPEED  OF  FASTEST  NEUTRONS  *  50  KM/MILLISEC. 

THIS  EFFECT  IMPORTANT  IN  MILLISECOND  TIME  FRAME. 

NO  THEORY.  CODE:  HAPS;  COULD  BE  INCLUDED  IN  LHAP. 

1.4  (HYDRODYNAMIC  MOTIONS)  x  (GEOMAGNETIC  FIELD) - ►EMF - ►CURRENTS  IN 

SURROUNDING  PLASMA — ►MAGNETIC  FIELD  PERTURBATIONS. 

THESE  EFFECTS  ARE  CALLED  MHD  EMP. 

IMPORTANT  FOR  LONG-LINE  GROUND  SYSTEMS. 

TIME  FRAME:  TENS  TO  HUNDREDS  OF  SECONDS. 


JSECKDU0  HM  BUMMiOV  IUMED 
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SEQUENCE  OF  EVENTS  RELEVANT  TO  MHD  EMP 
T  *  FEW  MILLISEC 


VdC/-  BURST 


•  GAMMA  RAYS  ABSORBED  BETWEEN  20  AND  50  KM. 

MAKE  EMP,  DURATION  ~1  y  SEC.  INDUCED  CONDUCTIVITY  DIES  AWAY  IN  0.1 
TO  1  MILLISEC. 

•  X  RAYS  ABSORBED  BETWEEN  70  AND  100  KM. 

MAKE  PLASMA  OSCILLATIONS.  INDUCED  CONDUCTIVITY  LINGERS  FOR  »  10  SEC, 
SCREENS  GROUND  OBSERVERS  FROM  EARLY  PART  OF  MHD  EMP. 


»»  ti*m**‘  ■+*****?*•  -»,<r  »  -».n*~  i'T)rvi>i  w 


2.2  T  *  0.1  SEC 


•  EXPANDING  BOMB  DEBRIS  AND  SHOCK  WAVE  PUSH  GEOMAGNETIC  FIELD  AWAY 
FROM  8URST  REGION. 


•  GEOMAGNETIC  PERTURBATIONS  SCREENED  FROM  GROUND  OBSERVERS  BY  CONDUCTING 
X-RAY  PATCH. 


•  SHOCK  WAVE  HEATS  A  LARGE  VOLUME  OF  THE  UPPER  ATMOSPHERE. 

•  8-RAY  PATCH  COULD  MAKE  MAGNETIC  PERTURBATIONS. 

PROBABLY  SMALL;  NO  EXPERIMENTAL  DATA. 
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2.3  T  »  1  SEC 


400  KM 


300  KM 


200  KM 


100  KM 


0 


GROUND 


•  BOMB  DEBRIS  AND  SHOCK-HEATED  AIR  AT  HIGH  ALTITUDES  MOVE  UPWARDS  AND 
OUTWAROS,  CARRYING  GEOMAGNETIC  FIELD  ALONG,  CAUSING  VERY  LARGE 
GEOMAGNETIC  FIELD  DISTORTIONS. 


•  THESE  PERTURBATIONS  STILL  SCREENED  FROM  GROUND  OBSERVERS  BY  CONDUCTIVITY 
IN  X-RAY  PATCH. 


•  SHOCK-HEATED  AIR  IONS  AT  LOWER  ALTITUDES  SPIRAL  DOWN  FIELD  LINES, 
ARE  STOPPED  IN  AIR  BETWEEN  100  AND  150  KM,  CREATE  KINETIC  ENERGY 
PATCH  (SOMETIMES  CALLED  DEBRIS  PATCH). 
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•  CONDUCTIVITY  IN  X-RAY  PATCH  IS  GONE  BY  RECOMBINATION. 


•  BOTTOM  OF  KINETIC  ENERGY  PATCH  EXPANDS  DOWNWARD,  COMPRESSING 
HORIZONTAL  COMPONENT  OF  GEOMAGNETIC  FIELD  AGAINST  EARTH. 


TOP  OF  KINETIC  ENERGY  PATCH  EXPANDS  UPWARDS  BUT  THIS  EFFECT  IS 
SCREENED  FROM  GROUND  OBSERVERS  BY  CONDUCTIVITY  IN  PATCH. 


GROUND 


ENTIRE  KINETIC  ENERGY  PATCH  EXPANDS  AND  RISES,  TENDING  TO  CARRY 
GEOMAGNETIC  FIELD  WITH  IT. 


HORIZONTAL  COMPONENT  OF  FIELD  BETWEEN  PATCH  AND  GROUND  IS  EXPANDED 
AND  REDUCED. 


TIME  (SECONDS) 

MAGNETOMETER  DATA  FROM  JOHNSTON  ISLAND 
FOR  FISHBOWL  DATA  TESTS 

1  GAMMA  =  10"5  GAUSS 


ELECTRIC  FIELD  PATTERNS 


I  *  m.  •***£ 


5.  FIELD  MAGNITUDES  EXPECTED 

•  STARFISH  6B  AT  JOHNSTON  ISLAND  WAS  -700  GAMMAS. 

FIELD  CHANGE  WAS  PROBABLY  CONSIDERABLY  LARGER  NORTH  OF  J.I. 

•  700  GAMMAS - *-E  «0.05  VOLTS/ KM  IN  OCEAN 

E  »  3.  VOLTS/ KM  IN  LAND 

•  MRC  HAS  SUGGESTED  THAT  E  MAY  BE  AS  LARGE  AS  30  VOLTS/ KM  FOR 
SINGLE  BURST,  100  VOLTS/KM  FOR  MULTIPLE  BURSTS  WITHIN  100 
SECOND  INTERVAL. 

•  BASIS  FOR  THESE  ESTIMATES  IS  EXCEEDINGLY  SHAKY  AT  PRESENT. 

•  SIZE  OF  GROUND  AREAS  EXPOSED  IS  SEVERAL  HUNDRED  KILOMETERS. 
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6.  VULNERABILITY  MECHANISMS 


•  PRIMARILY  OF  CONCERN  FOR  LONG-LINE  SYSTEMS. 


•  COPPER  WIRE  WITH  1  CM2  CROSS  SECTION,  TOO  KM  LONG  HAS  R  *  1 7  OHMS, 
SEES  VOLTAGE  3000  VOLTS.  IF  GROUNDED  AT  BOTH  ENDS.  CURRENT  ~  180 
AMPS  FOR  100  SECONDS.  HEAT  DEVELOPED  <s  5  JOULE/ CM^;  TEMPERATURE 
RISE  «  1 . 5°C.  NOT  BAD. 


•  IF  C0MPARA8LE  TOTAL  ENERGY  (107  JOULES)  WERE  DUMPED  INTO  A  PROTECTIVE 
DEVICE  AT  ONE  END  OF  THE  WIRE,  IT  COULD  BE  SERIOUS. 


•  IN  TYPICAL  LONG-LINE  COMMUNICATIONS  SYSTEMS  THE  WIRE  IS  NOT  GROUNDED 
(TO  PREVENT  TELLURIC  CURRENTS).  THEN  3000  VOLTS  APPEARS  BETWEEN 
EQUIPMENT  AND  GROUND.  MANY  SYSTEMS  THEN  TURN  THEMSELVES  OFF.  IN 
SOME  SYSTEMS,  D.C.  CURRENT  THROUGH  WIRE,  WITH  GROUND  RETURN,  IS 
USED  TO  POWER  REPEATERS.  THESE  MUST  BE  TURNED  OFF.  MAGNETIC  STORMS 
COMMONLY  CAUSE  SUCH  OUTAGES. 
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